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The Lya line of atomic hydrogen is often a luminous 
component of the radiation emitted by distant galaxies. Ex- 
cept for those galaxies which have a substantial central 
source of non-stellar ionizing radiation, most of the Lya 
radiation emitted by galaxies is generated within regions of 
the interstellar medium which are photoionized by starlight. 
Conversely, much of the energy radiated by photoionized 
regions is carried by the Lya line. Only hot, massive stars 
are capable of ionizing hydrogen in the interstellar medium 
which surrounds them, and because such stars are neces- 
sarily short-lived, Lya emission traces regions of active star 
formation. 

Observations of Lya radiation - including indirect ob- 
servations of thi effects of such radiation upon the gas in 
which it is generated - may provide a valuable probe of 
the interstellar medium in star-forming regions of external 
galaxies; in the case of very distant galaxies Lya observa- 
tions may be almost the only such probe. Careful theoretical 
analysis is required to interpret these observations, because 
the transport of H Lya photons is profoundly affected by 
the fact that in traversing an astrophysical medium such 
photons may suffer frequent absorption and resonant re- 
emission by hydrogen atoms in the ground state. Indeed, in 
escaping a neutral interstellar cloud of typical dimensions, 
H Lya photons may suffer millions of resonant scatterings. 
The transfer of Lya radiation in a highly opaque medium is 
thus a random walk process in which photons diffuse both 
in space and - by virtue of, the velocity dispersion of the 
scattering H atoms - in frequency. Diffusion in frequency 
results in a broadening of the line profile: in traversing a 
medium of atomic column density N(H°), Lya photons ac- 
quire a linewidth of 227 (10 -20 cm 2 .A(i7°)) 1 / 3 At;5^ 3 km s -1 
FWHM, where Av$ km s _1 is the velocity dispersion of the 
scattering H atoms. 

We argue that the strength of the Lya emission ob- 
served from external galaxies may be used to estimate quan- 
titatively the dust content of the emitting region, while the 
Lya line profile is sensitive to the presence of shock waves. 
Interstellar dust particles and shock waves are intimately 
associated with the process of star formation in two senses. 
First, both dust particles and shock waves owe their ex- 
istence to stellar activity; second, they may both [Serve} as 
agents which facilitate the formation of stars, shocks by trig- 
gering gravitational instabilities in the interstellar gas that 
they compress, and dust by shielding star-forming molecular 
clouds from the ionizing and dissociative effects of external 
TJV radiation. By using Lya observations as a probe of the 
dust content in diffuse gas at high redshift, we might hope 
to learn about the earliest epochs of star formation. 

Furthermore, the effects of Lya radiation in selectively 
pumping warm H 2 molecules may be inferred from the in- 


frared spectrum of at least one external galaxy; such effects 
may serve as signpost of interstellar gas which has been 
heated and partially dissociated and ionized by X-rays. 

II. Lya as a probe of interstellar dust 

In addition to being scattered, Lya photons may be de- 
stroyed by dust absorption, a process which dominates the 
heating of dust in HII regions. Indeed, because of the extra 
path length that they must travel, Lya photons are particu- 
larly vulnerable to attenuation by interstellar dust. The ef- 
fects of such attenuation have been described quantitatively 
by numerical solution of the transfer equation (Hummer and 
Kunasz 1980), and recently (Neufeld 1989) by an extension 
of an analytic method due to Harrington (1973) in which 
the essential features of the radiative transfer process are 
described by a Poisson equation. Figure 1 shows the frac- 
tion of Lya photons, f c , which can escape from the center 
of a dusty interstellar slab. 



The horizontal axis gives the atomic column density 
N(H°) from the center of the cloud to its surface, 
and the different curves are labelled with the quantity 
logioKdAvg 2 ^^ 0 )- 1 ], where x(H°) is the neutral atomic 
fraction, and fa is the dust opacity per hydrogen nucleus, 
normalized relative to the standard Galactic value of Draine 
and Lee (1984). 

In using observations of Lya line strengths to estimate 
the dust content in an external galaxy, we require an esti- 
mate of the intrinsic, unattenuated luminosity. If the in- 
trinsic equivalent width is known (from modelling of the 
stellar population), such an estimate may be derived from 
continuum measurements close to the Lya wavelength. Al- 
ternatively, observations of Balmer or Brackett lines, or of 
free- free radio emission, might be used to determine the 
emission measure of the ionized gas. The distribution of 
ionized gas and of Lya emission is also crucial. 

Hartmann et al. (1988) have recently presented mea- 
surements of the Lya/H/? ratio in several Markarian galaxies 
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of redshift ,js <; O.P."TBe results are plotted in Figure 2 as a 
function of the estimated [0]/[H] ratio. 



Without attenuation, the Lya/H/3 ratio is expected from 
the theory of radiative recombination to lie around 20; the 
observed values range from about 1 (the detection limit in 
this data set) to about 10, presumably reflecting the fact 
that the resonantly scattered Lya photons are far more vul- 
nerable to dust absorption than the 11/3 photons. Further- 
more, the ratio tends to be larger in galaxies of low metal- 
licity in which the dust content might be expected to be 
smaller. 

We have computed the expected Lya/H/3 ratio for the 
idealized case of a plane-parallel slab in which the neutral 
fraction, dust content, and recombination line emissivity are 
constant. The results are shown in Figure 3 for slabs of var- 
ious surface- to-center column densities N(H°). The scat- 
tering H atoms are assumed to have a velocity dispersion 
of 10 km s -1 , and each curve is labelled with the value of 
logi 0 (/V(i/ 0 )/cm- 2 ). 



The ratio decreases with increasing N(H°) until the slab 
becomes optically thick for absorption both of Lya photons 
and of H/3 photons. Further work is needed to treat the 
case where recombination line photons are generated within 
regions of low neutral fraction (HII regions) that are em- 
bedded in a medium of large neutral fraction (an interstellar 
cloud). 

Measurements of the Lya line width might in princi- 
ple provide a further constraint upon the dust content of an 
emitting galaxy, since the line width yields an upper limit on 
the atomic column density through which a typical Lya pho- 


ton has travelled. For example, Hunstead and Pettini (1989) 
have recently reported the detection of a narrow emission 
line from the hydrogen cloud which gives rise to the z=2.465 
damped Lya line in the spectrum of the QSO 0836+113. 
The narrowness of the line (< 50 km s -1 FWHM) implies 
that the detected Lya radiation has travelled through at 
most an atomic column of (1 x 10 18 /A«5) cm -2 . However, 
the total column density through the cloud is known (from 
the absorption line it produces) to be 4.2 x 10 2O cm -2 . Thus 
the Lya emission is either generated close to the near surface 
of the cloud and nowhere else, or we are prevented from see- 
ing embedded sources of radiation due to the effects of dust 
attenuation. A detailed quantitative analysis needs still to 
be carried out. 

III. Lya as a probe of interstellar shocks 

In considering the transfer of resonance-line radiation 
in a shocked interstellar cloud, we have found (Neufeld and 
McKee 1988) that the repeated scattering of Lya radia- 
tion across a shock front results in a systematic blueshift 
which may greatly exceed the shock velocity. The blueshift- 
ing process is a first-order Fermi process analogous to that 
invoked to explain the acceleration of cosmic rays. When 
Lya radiation is generated within a shocked cloud of neu- 
tral material the emergent line profile is given by: oc 

u 2 7 exp(— 0.034 u® 7 /uj7A 2 o), where = 100uj7 km s -1 is 
blueshift relative to line center, 100u S 7 km s -1 is the shock 
velocity, and N = ltP’A^o cm -2 is the column density of 
hydrogen atoms on either side of the shock front. The spec- 
trum peaks at a blueshift 270(Af 2O u J 7) 1 ^ 3 km s _1 . 

Fermi acceleration may be responsible for the blue 
asymmetric line profiles that have been observed in high 
redshift Lya galaxies. The Lya line profile in 3C326.1, mea- 
sured by Djorgovski (1988), is shown in Figure 4. 



Velocity shift (km/s) 


The observed spectrum may be accounted for quantita- 
tively by a model in which shocks, driven into a population 
of interstellar clouds by a radio lobe, trigger the formation of 
ionizing stars and Fermi accelerate the Lya radiation emit- 
ted by HII regions surrounding those stars. To match the 
observed line profile, which shows a blue wing extending to 
2000 km s -1 from line center, the required pressure behind 
the shocks is 6 x 10 -lo dyn cm -2 , consistent with a lower 
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limit on the pressure in the radio lobe that is set by the 
observed synchrotron emissivity. The shock velocity must 
lie between 10 and 400 kms' 1 and the pre-shock density 
between 0.16 and 260 cm -3 . The virial temperature of the 
clouds is 2 x 10 4 K, implying a line-of-sight velocity dis- 
persion of at least 17 km s -1 for thermal or turbulent sup- 
port. To allow the escape of Lya radiation, the dust opacity 
per hydrogen atom must lie at least a factor 60 below the 
standard Galactic value. The absolute strength of the ob- 
served Lya emission requires that the shocks trigger star 
formation with an efficiency greater than 0.5%. 

A further consequence of the Fermi acceleration mech- 
anism is that galaxy mergers, particularly between galaxies 
with low dust content, should produce Lya lines with strong 
blue wings. 

IV. H 2 pumping and the spectrum of NGC 6240 

In warm gas containing H 2 in vibrationally-excited 
states, H Lya radiation may pump coincident lines of the H 2 
Lyman bands (Shull 1978, Black and van Dishoeck 1987), 
thereby being converted into ultraviolet fluorescent photons 
with a characteristic spectrum that has been observed both 
in sunspots (Jordan et al. 1977) and in Herbig-Haro ob- 
jects (Brown et al. 1981; Schwartz 1983). There are two 
H 2 Lyman band transitions that lie close to the Lya rest 
frequency: the B-X 1-2 R(6) transition (at a redshift of 
15 km s _1 from line center) and the B-X 1-2 P(5) transition 
(at a redshift of 99 km s -1 .) The relative pumping rates in 
the two transitions may serve, in principle, as a probe of the 
Lya line profile in the region of the warm H 2 . 

The pumping process may be treated analytically by 
methods described by Neufeld (1989). Figure 5 shows the 
fraction of Lya photons, (1 — /„), which are converted into 
fluorescent photons of the Lyman bands before escaping 
from the center of a slab with a molecular fraction rela- 
tive to H atoms of x(H?) = n(// 2 )/n(//°). Each curve is 
labelled with the quantity — log 10 z:(J72)' 



which the H 2 energy levels are thermally populated and the 
H atoms and H 2 molecules have a thermal velocity disper- 
sion. Comparing Figures 1 and 5, we find that under these 
conditions, H 2 absorption dominates dust absorption when 
the molecular fraction exceeds about 10 -4 . Figure 6 com- 
pares the total rate of pumping in the two transitions. The 
vertical axis shows the number of pumps in the B-X 1-2 P(5) 


transition divided by the number in the B-X 1-2 R(6) tran- 
sition. 



log 10 [n(H z )/n(H°)] 


In addition to generating UV fluorescence, the pumping pro- 
cess may also modify the infrared emission spectrum of the 
warm, pumped H 2 . For example, we expect the (v=2, J=5) 
and (v=2, J=6) states o f H 2 to be selectively depopulated in 
regions of high Lya flux. Such depopulation is the probable 
explanation of the anomalously small v=2-l S(3) line inten- 
sity reported recently by Lester, Carr and Harvey (1988) for 
the starburst galaxy NGC 6240. (We might expect v=2- 
1 S(4) to be similarly diminished.) The required spatial 
coincidence of warm H 2 and Lya emission may argue in fa- 
vor of the model of Draine and Woods (1989), in which the 
copious H 2 rovibrational emission observed in NGC 6240 is 
produced in molecular clouds which are warmed, and par- 
tially ionized and dissociated by X-rays generated in a pop- 
ulation of embedded supernova remnants. The expected 
diminution in the v=2-l S(3) and v=2-l S(4) line strengths 
is currently being investigated quantitatively by Neufeld, 
Draine and Woods. 
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